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We analyse small-amplitude oscillations of a weakly viscous electrically conducting liquid 
drop in a strong uniform DC magnetic field. An asymptotic solution is obtained showing 
that the magnetic field does not affect the shape eigenmodes, which remain the spherical 
harmonics as in the non-magnetic case. Strong magnetic field, however, constrains the 
liquid flow associated with the oscillations and, thus, reduces the oscillation frequencies 
by increasing effective inertia of the liquid. In such a field, liquid oscillates in a two- 
dimensional (2D) way as solid columns aligned with the field. Two types of oscillations 
are possible: longitudinal and transversal to the field. Such oscillations are weakly damped 
by a strong magnetic field - the stronger the field, the weaker the damping, except for the 
axisymmetric transversal and inherently 2D modes. The former are overdamped because 
of being incompatible with the incompressibility constraint, whereas the latter are not 
affected at all because of being naturally invariant along the field. Since the magnetic 
damping for all other modes decreases inversely with the square of the field strength, vis- 
cous damping may become important in a sufficiently strong magnetic field. The viscous 
damping is found analytically by a simple energy dissipation approach which is shown 
for the longitudinal modes to be equivalent to a much more complicated eigenvalue per- 
turbation technique. This study provides a theoretical basis for the development of new 
measurement methods of surface tension, viscosity and the electrical conductivity of liq- 
uid metals using the oscillating drop technique in a strong superimposed DC magnetic 
field. 



1. Introduction 

Shape oscillations of levitate d metal droplets can be used to measur e the surface tension 



and viscosity of liquid metals (jRhim et al. I I1999L lEgry et al. 2005) . Theoretically, the 



former determines the frequency, while the latter accounts for the damping rate of oscil- 
lations. In the reality, experimental measurements may be affected by several side-effects. 
Firstly, levitated drops may be significantly aspherical and the oscillations amplitudes 
not necessarily s mall, whereas the clas sical theories describing the oscillation frequencies 
I Ravleigh 1945 ) and damping rates ( Lamb |[l993t IChandrasekhar 1981 ; Reid 1960h 
assume small-amplitude oscillations about an ideally spherica l equilibrium shape. Cor 



rectio ns du e to the drop asphericity have been c a lculated bylCummings fe Blackburn 



( 1991 ) and Survanaravana fe Bavazitoglu ( 199ll ) . iBratz fe Egrv ( 19951) find the same 
order correction to the damping rate resulting also from AC-magnetic field. The ef- 
fect of a moderate ampl i tude o n the oscillations of inviscid drops has been analysed by 
Tsamopoulos fe Brown ( 19831 ) who find that the oscillation fre quency decreases with 
the square of the amplitude. Using a boundary-integral method, Lundgren fc Mansour I 
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(|l99lh show that small viscosity has a relatively large effect on the resonant-mode cou- 
pling phenomena in the nonlinear oscillations of large axially symmetric drops in zero 
gravity. Nu merical simulatio n of large-amplitude axisymmetric oscillations of viscous liq- 
uid drop bv lBasaran (1992), who uses the Galerkin/finite-element technique, shows that 
a finite viscosity results in a much stronger mode coupling than predicted by the small- 
viscosity approximation. 

Secondly, the measurements may strongly be disturbed by AC-driven flow in the 
drop. The mode coupling by the internal circulation in axisym metrically oscillating drop 
has been studied numerically by iMashavek fc Ashgriz (1998) using the Galerkin/finite- 
element technique. To reduce the strength of the AC field necessary for the levitation 
and, thus, to minimise the flow, experiments may be conducted und er the microgravit y 
conditions during parabolic flights or on the board of space station l Egrv et al. I 19991 ). 
A cheaper alternative might be to apply a sufficiently strong DC magnetic field that can 
not only stabilise AC-driven flow but also suppress the convective heat and momentum 
transport resp onsible for the mode coupling under the terrestrial conditions as origi- 
nally shown by IShatrov et al. ( 20031) . Such an approach has been implemented first by 
Yasuda et al. I (|2004f ) on the electromagnetically levitated drops of Copper and Nickel 



which were submitted to a DC field of the induction up to 10 T. The only motion of Cu 
drops observed to persist in magnetic field above 1 T was a solid-body rotation about 
an axis parallel to the magnetic field. No shape oscillations, usually induced by the AC- 
driven flow fluctuations, were observed. Note that this implies only the suppression of 
AC-driven flow but not of the shape oscil l ations themselves which require an external 
excitation to be observable. I Yasuda et all ( 2005h study the effect of suppression of the 
melt flow on the structure of various alloys obtained by the electromagnetic levitation 
melting technique in a strong superimposed DC magnetic field. The use of high magnetic 
fields in various material processing applications is reviewed by lYasuda ] (l2007h . 

Note that a simila r suppression of AC- driven flow can also be achieved by a fast 
spinning of the drop (IShatrov et al. I [2 0071 that may be driven by an electromagnetic 
spin- up instability (jPriede fc Gerbeth II2000L 120061 ) . The effects of both the drop rotation 
and AC-driven flow on the frequency spectrum of shape osc i llatio ns have been modelled 
numerically bv lBoiarevics fc Pericleous l (l2009l ). IWatanabel(l2009l) demonstrates numer- 
ically that a large enough oscillation amplitude can compensate for the effect of rotation 
on the frequency shift. 

A novel method of measuring thermal conductivity of liquid silicon using the elec- 
tromagnetic levi t ation in a strong superimposed DC mag netic has been introdu ced by 
Kobatake et al. 1 2007 ). Subsequent numerical modelling bv lTsukada et al. 1 20091 ) shows 
that applying a DC magnetic field of 4T can suppress convection in molten silicon 
droplet enough to measure its real thermal conductivity. Later on this method has 
been extended to the m easurements of heat capacity of molten austeni tic stainless steel 
llFukuvama et al 2009) and also that of supercooled liquid silicon ( Kobatake et al. 



2010). 



In order to determine the surface tension and viscosity or the electrical conductivity 
one needs to relate the observed surface oscillations with the relevant thermophysical 
properties of the liquid. General small-amplitude shape oscil lations o f cond ucting drop 
in a uniform DC magnetic field have been analysed first by Gailit isl (jl966l ). Although 
a magnetic field of arbitrary strength is considered, the solution is restricted to inviscid 
drops. Moreover, only the frequency spectrum and magnetic damping rates are found 
but not the associated shape eigenmodes, which may be useful for experimental iden- 
tification of the oscillation modes. Energy dissipation by ax isymmetric oscil lations of a 
conducting drop in a weak DC magnetic field is considered bv lZambran (1966), who finds 
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Figure 1. Sketch to the formulation of problem. 



the magnetic damping rates in agreement with more general results of Gailitis ( 19661 ). 



Axisymmetric oscillations of an electromagnetically levitated drop of molten Al in a 



superi mposed DC magnetic field are modelled numerically by iBoiarevics fe Pericleous 



(j2003l ). A moderate DC magnetic held is shown to stabilise AC-driven flow and, thus, to 
eliminate the associated shape oscillations. A three-dimensional numerical simulation of 
an oscillating liqui d metal drop in a uniform static magnetic field has been carried out 
by iTagawa" ( 2007 ). The numerical results show that vertical magnetic field effectively 



damps the flow, while horizontal field tries to render the flow two-dimensional. 

In the present paper, we analyse free oscillations of a viscous e lectri cally conducting 
drop in a homogeneous DC magnetic field. In contrast to iGailitis (1966), we assume the 
viscosity to be small but non-zero and the magnetic field to be strong. This allows us 
to obtain an asymptotic solution to the eigenvalue problem for gener a l sma ll-amplitude 
3D shape oscillations including the eigenmodes left out by IGailitis I (|l966l) . which are 
necessary for the subsequent determination of the viscous damping. Firstly, we show that 
the eigenmodes of shape oscillations are not affected by strong magnetic field. Namely, 
they remain the spherical harmonics as in the non-magnetic case. The magnetic field, 
however, changes the internal flow associated with the surface oscillations and, thus, 
the frequency spectrum. As the drop oscillates in a strong magnetic field, the liquid 
moves as solid columns aligned with the field. Two types of such oscillations are possible: 
longitudinal and transversal to the magnetic field. The oscillations are weakly damped by 
a strong magnetic field, except for both the axisymmetric transversal and inherently 2D 
modes. The former are magnetically overdamped because the incompressibility constraint 
does not permit an axially uniform radial flow. The latter, which are transversal modes 
defined by the spherical harmonics with equal degree and order, I = m, are not affected 
at all because these modes are naturally invariant along the field. Because the magnetic 
damping for all other modes decreases inversely with the square of the field strength, the 
viscous damping may become important in a sufficiently strong magnetic field. 

The paper is organised as follows. The problem is formulated in <J21 Section [3] presents 
an inviscid asymptotic solution which yields the shape eigenmodes and frequency spec- 
trum of longitudinal and transversal oscillations. Magnetic damping is found in ^3.21 as a 
next-order asymptotic correction to the frequency. Viscous damping rates are calculated 
in ^3.31 first by the eigenvalue perturbation technique for the longitudinal modes and 
then by an energy dissipation approach for both of the oscillation modes. The paper is 
concluded by a summary and discussion of the results in S|4] 
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2. Problem formulation 

Consider a spherical non- magnetic drop of an incompressible liquid with radius Rq, 
density p, surface tension 7, electrical conductivity a, and a small dynamic viscosity 77 
performing small-amplitude shape oscillations in a strong uniform DC magnetic field B 
as illustrated in figure [TJ The velocity of liquid flow v and the pressure distribution p are 
governed by the Navier-Stokes equation with electromagnetic body force 

pd t v = -Vp + r]V 2 v + j x B, (2.1) 

where the induced current follows from Ohm's law for a moving medium 

j = a(E + v x B). (2.2) 

Owing to the smallness of oscillation amplitude, the nonlinear term in (|2.ip as well as 
the induced magnetic field are both negligible. In addition, the characteristic oscillation 
period tq is supposed to be much longer than the magnetic diffusion time /iquRq, where 
po is the permeability of vacuum. This leads to the quasi-stationary approximation ac- 
cording to which V x E = and E = — Vy, where ip is the electric potential. The 
incompressibility constraint V • v = and the charge conservation condition V • j = 
applied to (|2.ip and (|2.2p result, respectively, in 

V 2 p = a(B- V)(B-v), (2.3) 
V 2 tp = B-V xv. (2.4) 

For a uniform B under consideration here, applying the operators V x Vx, (B ■ V)B- 
and (B ■ V)B • Vx to (|2~T1) and taking into account V x V x (j x B) = a(B ■ V) 2 v 
together with (|2.3p and (|2.4p . we obtain 

[ P V 2 d t + a(B ■ V) 2 - 77V 4 ] {p, v} = 0. (2.5) 

Although the equation above applies to p, (p and v separately, these variables are not 
independent of each other. Firstly, owing to the incompressibility constraint, only two 
velocity components are mutually independent. Secondly, velocity is related to the pres- 
sure and electric potential by (|2.1j) , which can be used to represent v in terms of p and 
ip as done in the following. 

Boundary conditions are applied at the drop surface S defined by its spherical radius 
R = Ra + Ri(9, (f>, t), where R± is a small perturbation, which depends on the poloidal 
and azimuthal angles, 8 and (j), and the time t. The radial velocity at the surface is related 
to the radius perturbation by the kinematic constraint 

v R \ s = d t R 1 . (2.6) 

The normal component of the current at the drop surface, which is assumed to be sur- 
rounded by vacuum or insulating gas, vanishes, i.e., j n \ s = 0. In addition, there is no 
tangential stress at the free surface: 

n ■ d T v + t ■ d n v\ s = 0, (2.7) 

while the normal stress component is balanced by the capillary pressure 

Po + p - 2-qd n v n = 7V • n, (2.8) 

where po — 2j/Rq is the constant part of pressure, r is a unit tangent vector and n = 
V(i? — i?i)/|V(i? — i?i)| is the outward surface normal. For small-amplitude oscillations 
defined by R\ <§; Rq, we have n w — Vi?i. 

Henceforth, we proceed to dimensionless variables by choosing the radius Rq and the 
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characteristic capillary pressure Pq = •j/Rq as the length and pressure scales. The char- 
acteristic period of capillary oscillations is determined by the balance of inertia and 
pressure, which yields the time scale To = y R^ppy. The velocity and potential scales are 
chosen as vq — Rq/tq and <pa = vqBRq, respectively, where B = \B\ . In the dimensionless 
variables, (|2.5[) takes the form 

[V 2 d t + Cm(e ■ V) 2 - CaV 4 ] {p, tp, v} = 0, (2.9) 

where e = B/B is a unit vector in the direction of the magnetic field and Ca = i]/\/Rofry 
and Cm — aB 2 R^/ 'y/Ropj are the conventional and magnetic capillary numbers, respec- 
tively. They are the ratios of the capillary oscillation time tq defined above and the viscous 
and magnetic damping times, which are t v — pR^/v and r m = p/(o~B 2 ), respectively. In 
the dimensionless form, the normal stress balance condition (|2.8j) reads as 

(V 2 +2)R 1 +p-2Cad R v R \ R=1 = 0. (2.10) 

In the following, we assume viscosity to be small but the magnetic field strong so that 
Ca -C 1 and Cm S> 1, which means that the second and third terms in (|2.9p are much 
greater and much smaller, respectively, than the first one. Thus, we first focus on the 
effect of the magnetic field and ignore that of viscosity, which is considered later in £13.31 



3. Inviscid asymptotic solution 

Here we ignore viscosity that allows us to formulate the problem in terms of p, if and 
R\. Projecting the dimensionless counterpart of (|2.1I) . which takes the form 

Cmv + d t v = -Vp + CaV 2 v + Cm [ex Vp + e(e • «)] , (3.1) 

onto e R and e, and putting Ca = 0, we obtain 

CmvR + d t VR = —eR ■ Vp + Cm [e^ x e ■ + e# • ev n ] , (3.2) 
d t v u = -e ■ Vp, (3.3) 

where v„ = e • v is the velocity component along the magnetic field. Differentiating (|3.2H 
with respect to t and substituting d t v u from (|3.3p . we represent (|2.6p in terms of p and <p 

Cmd 2 t R x + d*Rx = [Cm (e R x e • Vd t f ~ (e R ■ e)e ■ Vp) - e R ■ Vd tP ] R=1 . (3.4) 

Velocity has to be eliminated also from the electric boundary condition given by the 
radial component of Ohm's law 

j R [ R=1 = -e R ■ [Vip + e x v] R=1 = 0. (3.5) 

Firstly, applying (Cm + dt) to (|3.5p and then using (|3.ip . we obtain 

[Cm(e R ■ e)e • Vip - e R x e • Vp + e R ■ Vd t (p] R=1 = 0. (3.6) 

In the inviscid approximation, (|2.I0p reduces to 

p|«=i = -(V 2 +2)i? 1 . (3.7) 

In the following, besides the spherical coordinates (R,9,<f>), we will be using also the 
cylindrical ones (r, cj>, z) with the axis aligned along the magnetic field so that e = e z . 

Solution is sought in the normal mode form {p, tp, Ri} = {p : <p, R\(r)e^ t+lm ^ , where j5, 
tp and R are axisymmetric amplitude distributions, m is the azimuthal wave number, and 
f3 is a generally complex temporal variation rate which has to be determined depending 
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on m, Cm and Ca. Then boundary conditions (|3.4|) . (|3.6[) and (J3T7J) for the oscillation 
amplitudes at R = 1 take the form 

f3 2 R + impcp + zd z p = -Cm- 1 (l3 3 R + (3d R p), (3.8) 
zd z (p = — Cm" (imp + ^8r0), (3.9) 
P=-(L Z + 2)R, (3.10) 

where L z = 4- ((1 — z 2 )-^) ~ IS the angular part of the Laplace operator in the 
spherical coordinates for the azimuthal mode m written in terms of z = cos 9. Further, 
it is important to note that 

L z pr(z) = -i(i+i)pr(z), (3.H) 

where P™(z), the associated Legend re function of degree I and ord er to, is an eigcnfunc- 



tion of L z with eigenvalue — 1(1 + 1) (jAbramowitz fc Stegun 1119721 ). Equation (|2.9p for p 
and <p can be written as 

[d 2 z + Cm-\L r + d 2 )((3 - Ca(L r + d 2 ))] {p, = 0, (3.12) 

where L r = d 2 +r~ 1 d r — m 2 /r 2 is the radial part of the Laplace operator in the cylindrical 
coordinates for the azimuthal mode m. Here we put Ca — 0, suppose Cm 3> 1, and search 
for an asymptotic solution in the terms of a small parameter Cm" 1 as 

{f>, <P,R,P} ~ {po,^o,-Ro, A)} + Cm~ l {pi,<pi,Ri,P\} H . 

Note that although (|3 . 1 2|) admits solutions with f3 ~ Cm found by iGailitis I (Il966l) . 



such quickly relaxing modes cannot be related with the surface deformations. From the 
physical point of view, drop is driven to its equilibrium shape by the surface tension, 
and the magnetic field can only oppose but not to accelerate the associated liquid flow. 
From the mathematical point of view, /3 ~ Cm ^> 1 applied to (|3 . 8[) results in Rq = 0, 
which means no surface deformation at the leading order in agreement with the previous 
physical arguments. Consequently, these fast modes represent internal flow perturbations 
which are not relevant for the shape deformations under consideration here. 

3.1. Oscillation frequencies 
At the leading order, (|3.12j) reduces to d 2 {jjq, tpo} = 0, whose general solution is 

{Po, fio}{r, z) = {p+, tp$}(r) + z{p~Q , <$>q }(r), (3.13) 

where the first pair of particular solutions are the functions of r only, while the second pair 
is linear in z but general in r. Owing to the z-reflection symmetry of the problem these 
two types particular of solutions do not mix and, thus, they are subsequently considered 
separately. We refer to these solutions in accordance to their z-parity as even and odd 
ones using the indices e and o. As shown below, the odd and even solutions describe 
longitudinal and transversal oscillation modes, respectively. 

3.1.1. Longitudinal modes 

For the odd solutions {pg, (fo}(r, z) = z{p^ , <^q }(r), boundary condition (|3.9|) . which at 
the leading order reads as zd z tpo = 0, results in <^q (r) = 0. The two remaining boundary 
conditions p. 81) and (|3.10p take the form 

(ifRZ = -zpo, (3.14) 
(L z + 2)R° = -zp^. (3.15) 
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Figure 2. Shapes and the associated liquid oscillations in the (a;,z)-plane parallel to the mag- 
netic field for the first four longitudinal oscillation modes with indices (I, m) = (2,1) (a), (3,0) 
(b), (3,2) (c), and (4,1) (d). 

Eliminating the pressure term between the equations above, we obtain an eigenvalue 
problem in /3q for Rq 



(L z + 2- (3° 2 )R° = 0, 



which is easily solved by using (|3.11[) as 



R° (z) 
/3 ° 



R° P[ 



(*), 



±i^(i-m+2), 



(3.16) 

(3.17) 
(3.18) 



where Rq is a small amplitude of oscillations and I — m is an odd positive number. Note 
that imaginary /3§ describes constant-amplitude harmonic oscillations with the circular 
frequency |/3 pl which diff ers from the corresponding non-magne tic result only by the 
factor of \fl (jLamb 1993 ). and coincides with the result stated by Gailitis ( 19661 ). Thus, 
strong magnetic field changes only the eigenfrequencies but not the eigenmodes of shape 
oscillations which, as without the magnetic field, are represent ed by separate spherical 



funct ions (associated Legendre functions with integer indices) (jAbramowitz h Stegun 
1972). Similarly to the non-magnetic case, the frequency spectrum for odd modes is 



degenerate because it depends only on the degree I but not on the order m. Thus, for 
each I, there are [1/2] odd modes with different m. 

Taking into account that z = \/\ — r 2 at the surface, the radial pressure distribution 
is obtained from (|3.14p as 



Po (r) = -^R^Vl-r^/VT^ 2 



(3.19) 



According to (|3.3[) , this pressure distribution is associated with the axial velocity com- 
ponent 

») = -/^Po M, (3.20) 
while two other velocity components transversal to the magnetic field are absent in the 
leading-order approximation. Thus, the liquid effectively oscillates in solid columns along 
the magnetic field as illustrated in figure[2]for the first four longitudinal oscillation modes 
defined by the indices (l,m) = (2,1), (3,0), (3, 2), and (4,1). Since such a flow does not 
cross the flux lines, the oscillations are not damped by the magnetic field in the leading- 
order approximation. 

3.1.2. Transversal modes 

For the even solutions {Pq, <^ol( r J z ) — {Po > < /?o~}( r )j the leading-order boundary condi- 
tion (|3.9p is satisfied automatically. The two remaining conditions (|3.8p and (j3.10p then 
take the form 



■ irrupt 



0, 



(3.21) 
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{L z + 2)R e = -p+. (3.22) 

In contrast to the longitudinal modes considered above, now we have two equations (|3.2ip 
and ()3.22j) but three unknowns. To solve this problem, we need to consider the first-order 
solution to (|3.12l) which now takes the form d 2 {p\ 1 <£f } = — /3oL r {pJ, (p$} and yields 

{p{, #}(r, z ) = {Pt> vtHr) ~ \^z 2 L r {p+, (3-23) 

Then boundary condition (13.91) results in imp^ — /3Q(z 2 L r — rd r )<p~Q = 0. Combining this 
with (|3.21[) and (|3.22l) and taking into account 

z 2 L r -rd r \ R=l = L z +m 2 , (3.24) 

we obtain 

[L z + 2 + (/3 e a /m) 2 (L z + to 2 )] R e = 0. (3.25) 
Using (|3.1ip . we readily obtain 

Rl{z) = R%Pr(zh (3.26) 

where Rq is a small oscillation amplitude and I — m is an even non- negative number . 
The result above again agrees with the asymptotic solution given by Gailitis ( 19661 ). 



Similarly to the odd solutions found in the previous section, even eigenmodes are repre- 
sented by separate spherical functions, and the oscillations are not damped at the leading 
order. In contrast to the odd modes as well as to the non-magnetic case, the frequency 
spectrum (|3.27[) is no longer degenerate and frequencies vary with the azimuthal wave 
number to. In particular, there are two important results implied by Q3.27p . Firstly, the 
oscillation frequency for the axisymmetric modes specified by to = is zero. This means 
that these modes are over-damped and do not oscillate at all. Secondly, the oscillation 
frequency for the modes with m = I is exactly the same as without the magnetic field, i.e., 
— + 2). This is because the liquid flow associated with these o scillation mode s 



is inherently invariant along the field and, thus, not affected by the last (jGailitis 11196 



The electric potential and pressure distributions follow from (|3.21|) and (|3.22ll as 

0§(r) = im-^i^Vl-r 2 ), (3.28) 



Po 



(r) = (l-l)(l + 2)R%(Vl-r 2 ). (3.29) 



The associated velocity distribution is obtained from (I3.1[) . Firstly, equation (|3.3p implies 
that the liquid oscillations are purely transversal to the magnetic field. In the leading- 
order terms, we obtain from (|3.f p 

^(r^) = e z xV#,ft (3.30) 

which shows that the velocity is not only transversal but also invariant along the magnetic 
field. Thus, the liquid again oscillates as solid columns, but in this case transversely to 
the field which has no effect on such a flow. This is because the e.m.f induced by the flow, 
which is invariant along the magnetic field, is irrotational, i.e., Vx (vxB) = (-B-V)u = 0, 
and, thus, unable to drive current circulation in a closed liquid volume. 

Note that for the axisymmetric modes (to = 0), the potential f|3 . 28[) and the associated 
velocity (|3 . 30(1 take an indeterminate form. Namely, for m = 0, boundary condition (|3.2ip . 
which in this is case straightforwardly implies a zero frequency, is satisfied by an arbitrary 
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Figure 3. Shapes and the associated liquid flows in the horizontal mid-plane (z = 0) perpen- 
dicular to the magnetic field for the first four transversal oscillation modes defined by indices 
(J, m) = (2, 2) (a), (3, 1) (b), (3,3) (c) and (4, 2) (d). 



potential distribution independent of the radius perturbation. As seen from p.28|) . a non- 
zero axisymmetric potential is associated with a purely azimuthal velocity. Consequently, 
this mode is irrelevant and can subsequently be neglected because it represents an internal 
flow perturbation which is just compatible but not coupled with axisymmetric shape 
deformations similarly to the fast modes discussed at the end of Moreover, this is 
consistent with (|3.27[) according to which axisymmetric transversal modes are static in 
the leading-order approximation that implies a zero velocity and, consequently, a zero 
associated potential. 

Expression (|3.30|) implies that the velocity streamlines coincide with the isolines of 
ipo, which, thus, represents a stream function for the flow oscillations. Figure [3] shows 
the shapes and streamlines of the associated liquid flow in the horizontal mid-plane for 
the first four transversal oscillation modes. Note that the first and the third mode with 
the indices (I, m) = (2,2) and (3,3) are both naturally invariant in the direction of the 
magnetic field and, thus, effectively non-magnetic. The second mode with (I, m) — (3, 1) 
corresponds to the drop oscillating in such a way that horizontal cross-sections remain 
circular in the small-amplitude limit under consideration while the whole shape deforms 
because of vertical offset of the cross-sections. 



3.2. Magnetic damping 

3.2.1. Longitudinal modes 

In order to determine the magnetic damping rates for longitudinal modes, we have to 
consider the first-order solution governed by 

d 2 z {p° 1 ,<p° 1 } = -P° zL r {pn,0}, 

which yields 

{p?,#}(r,z) = z{pi,0i}(r) - i/3 °z 3 L r {p -,0}. 

Then (|3.9[) and (I3.14[) applied consecutively result in z(pi — —imzpg, which combined 
with (EU, (jBTTUj) and (jXTlf yields 

(L z + 2 - pf)Rl = ±(3fz [z 2 L r - 3rd r ] z^R" + /3 °(2/?° - m 2 (3f)R°. 

After some algebra, we obtain z [z 2 L r — 3rd r ] z' 1 ]^^ = L z + 2 + m 2 , and, consequently, 



(L z + 2 - pf)Ri = -P° [pf {L z - 2m 2 + 2) + 6/3? 



(3.31) 
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The l.h.s. operator above is the same as that in ()3.16|) which has Rq as its eigensolution 
with a zero eigenvalue. Owing to (|3.11ll and (|3.17p . Rq is an eigensolution of the r.h.s 
operator of (|3.3ip . too. Thus, for (|3.3ip to be solvable, its r.h.s has to be free of the terms 
proportional to Rq, that yields 

PI = ~(l - l)(i + 2)((i - + 2) + 2m 2 ). (3.32) 

Note that conversely to the frequency for longitudinal oscillation modes (|3 . 1 8|) , the mag- 
netic damping rate above is not degenerate and varies with m. 

3.2.2. Transversal modes 

Similarly to the oscillation frequency considered above, boundary conditions (|3.10[) 
and p. 81) applied to the first-order solution (|3.23[) result in 



1 



/3 Rt + im<pf= pf \rn- 2 (z 2 L r -rd r )' z - ^z z L r -lj - PI Rq, (3.33) 

{L z + 2)R\ = -p\. (3.34) 
To solve this first-order problem, we again need a second-order solution governed by 

d 2 z {P2, <f>%) = ~P e (L r + - P1L r {p+, <p+}, 

which, by taking into account (|3.23[) . yields 



ftfaz) = 0+(r) ~ ^ 2 L r ^+ + p* 



4! 



Li 



Cp+ - -P{z 2 L r ^+, (3.35) 



Then (f3T0|) results in 

impl — /?g(z 2 .L r — rd r )ip^ = im~ 1 



P e Q{z 2 L r - rd r ) - -Pfz 2 (z 2 L r - ird r )L 7 
6 



Substituting (p\ and p\ from p.33p and (|3 . 36[) into (|3.34[) and using 

z 2 (z 2 L r - 3rd r )L r - 3(z 2 L r - rd r )z 2 L r \ R=1 = 2m 2 - 2(L Z 



Rq. 

(3.36) 



after some algebra we obtain an equation for Rf, which is the same as (|3.25[) for Rq, 
except for the r.h.s. that now reads as 



3_ 

3to 2 



[(P e /m) 2 ((L z + m 2 ) 2 - m 2 ) (3L Z + 2m 2 ) - 6pt(L z + m 2 )] R e 



By the same arguments as for (I3.31[) . the solvability condition applied to the expression 
above results in 



Pt = - 



(l - + 2)(l 2 - m 2 ){{l + l) 2 - m 2 ){3l(l + 1) - 2m 2 ) 



6(1(1 + 1) -m 2 ) 2 

which again coincides with the corresponding result of iGailitis" 1 1966I ). 



(3.37) 



3.3. Weak viscous damping 

There are three effects due to viscosity in this problem. Firstly, viscosity appears in 
the normal stress balance condition (|2.10[) as a O(Ca) correction to the inviscid solu- 
tion obtained above. Secondly, viscosity also appears as a small parameter Ca in (|3. 12|) 
which again implies the same order correction when the leading-order inviscid solution is 
substituted into this term. Thirdly, viscosity enters the problem implicitly through the 
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free-slip boundary condition (|2.7|) which was ignored by the inviscid solution but needs 
to be satisfied when viscosity is taken into account. To satisfy this condition, the leading- 
order solution needs to be corrected by the viscous term in (|3.12p . where Ca appears as 
a small parameter at the higher-order derivative. For this small viscous term to become 
comparable with the dominating magnetic term at the surface, the expected correction 
has to vary over the characteristic length scale 8 ~ y ' Ca/ Cm = Ha -1 , which is defined 
by the Hartmann number Ha = BqRq^Jct / (pv). Moreover, for the viscous correction of 
the tangential velocity v T in the Hartmann layer to compensate for a O(l) tangential 
stress due to the leading-order inviscid solution, v T ~ Ha^ 1 is required. Then the incom- 
pressibility constraint implies an associated normal velocity component of an order in 8 
smaller than v T , i.e., v n ~ Ha~ 2 . This normal velocity correction is subsequently neg- 
ligible. But this not the case for the tangential velocity correction v T , which according 
to (|2.3p is expected to produce a pressure correction p ~ Cm/ Ha 2 ~ Ca. The last is 
comparable with the normal viscous stress produced by the leading-order inviscid flow. 
Taking into account the estimates above and (p ~ Sv^ ~ Ha~ 2 , which follows from (j2.4[) . 
we search for a viscous correction as 

{p, tp,v} ~ {p 0l ip 0l v } + Ca{poi,ipoi,v i} + {Cap, Ha~ 2 tp, Ha~ l v] 
{R, [3} ~ {Ro, A>} + Ca{R 01 ,(3 01 } + • • ■ , 

where the terms with the tilde account for a Hartmann layer solution localised at the 
surface. 

3.4. Eigenvalue perturbation for longitudinal modes 

We start with the core region, where the additive boundary layer corrections are supposed 
to vanish. The first-order viscous corrections for the pressure and potential {Pqi, <^oi}( r J z) — 
z{pQi,fiQi}(r) are obtained similarly to the leading-order inviscid solution (|3 . 1 3|) . Now, 
instead of the kinematic and electric boundary conditions (|3.4jl and (|3.6| derived in the 
inviscid approximation, we have to use the original ones (|2.6j) and (|3 . 5[) containing the 
velocity, which again follows from the Navier-Stokes equation (|3. lj) including the viscous 
term ~ Ca. 

For the longitudinal modes, described by the odd solutions, (|3.1|) yields 

«i + /W = -Poi + L r w° , (3.38) 
u° ol = e z xD0° ou (3.39) 

where w and u are the velocity components parallel and perpendicular, respectively, to 
the field direction e z , and D = e - lm< A Ve lm(?i is a spectral counterpart of the nabla oper- 
ator for the azimuthal mode m. Since, as shown above, both the potential and velocity 
perturbations in the Hartmann layer are higher-order small quantities and, thus, negli- 
gible with respect to the core perturbation, the electric boundary condition ()3.5j) can be 
applied at R — 1 directly to the first-order core solution as dRip^ — ■ e z x u^. 
Taking into account (|3.39p . this yields 0q X = and, hence, tt^ = 0. Consequently, the 
first-order velocity perturbation in the core for the odd modes is again purely longitu- 
dinal. Then, the kinematic constraint (|2.6|) for the leading- and first-order terms takes, 
respectively, the form 

zw° = «, 
zw° 01 =(3° R° 01 +p° 01 R° . 
These expressions combined with (|3.38[) result in 

P° (P° R° 01 + 2/M - zL r z- l IQ = -p° 01 , (3.40) 



12 



J. Priede 



which defines the first-order core pressure perturbation at R = 1. In addition, we need 
also the Hartmann layer pressure correction which according to the estimates above is of 
the same order of magnitude as the core one. 

To resolve the Hartmann layer, we introduce a stretched coordinate R = (1 — R)/S 
l Hinch|[l99lh . where 8 = Ha 1 is the characteristic Hartmann layer thickness. In the 
Hartmann layer variables, (|3.12[) takes the form 

(Cm" 1 ft + z 2 - d\)d\{p, 0, i} = 0. (3.41) 

For Cm 3> 1, the inertial term ~ Cm^ 1 is negligible in (|3.4ip with respect to the magnetic 
one ~ z 2 , except for \z\ < Cm^ 1 ^ 2 . First, ignoring this term, which, as shown below, 
gives a next-order small correction, the solution of (|3.4ip vanishing outside the Hartmann 
layer can be written as 

fafrv} = v'}(z)e-M* (3.42) 

where the index s denotes the surface distribution of the corresponding quantity. Then 
the free-slip boundary condition (|2.7|) results in 



v% = -W^imr-^R + d R (v ^/R)), (3.43) 

v 6 e = -\z\-\devoM + d R (vo,o/R)). (3.44) 

For the longitudinal modes, defined by the odd solutions, the leading-order inviscid ve- 
locity is purely axial 

v° = e z w (r) = -e z !3° - 1 Po(r). (3.45) 
Substituting this into (|3.44[) and taking into account that the radial pressure distribution 
at the surface is related to the radius perturbation by (|3.14p , we obtain 

12 _ 2\ d fro 

v$=p K ' 3 ' 46 

z\z\ az z 

Pressure is related to the velocity by (|2.3p . which in the dimensionless form reads as 
V 2 p = Cmd z v z . In the Hartmann layer variables, this equation takes the form 

d\p = rzd k v e . (3.47) 

Substituting the general solutions for pressure and velocity given by ()3.42j) into (|3.47| 
and using ()3.46j) . we find 

p s = -rzlzl" 1 ^. (3.48) 
Substituting the normal component of viscous stress 

together with the core and boundary layer pressure contributions defined by (|3.40l) and 
(I3.48P into the normal stress balance condition (|2.10p . we finally obtain 



(L z + 2 - P° 2 )R° Q1 = /3° 



2p° m R° - z- 2 {L z + m 2 + 2)R° - 2(1 - ^ 2 )^^ 



(3.49) 



The sought for viscous damping rate is obtained in the usual way by applying the solv- 
ability condition to (|3.49[) that after some algebra results in 



-(Z(Z + l)-m 2 -2)^-^ (3.50) 
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m = 1 2 3 4 5 6 
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35 
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51 
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27 




154 
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44 


2 


22 




169 




403 




135 


2 


125 


fi 


95 



Table 1. The viscous damping rates — /?oi f° r the first 6 longitudinal oscillation modes. 



if (z) 



= -(21 + 1) 



2^+1) 



to 2 ) - 1 - If 



where 



Tin 

1 l — 



{I - m)\ 
{I -m)! 



^) , d jf co 



TO 



J£ a + (2Z 



dz 



dz 



1) 



p _ m 2 

can be calculated from the above recurrence relation starting with I 
into account that I\ a = for I < to. For the modes with to = I - 



(3.51) 

= TO + 1 and taking 
1, we have /3gi = 



^(2Z + 1)(7 — 1) , which is th e half of the corresponding viscous damping rate without the 
magnetic field (jLamb 1119931) . Although the viscous damping rate increases for smaller to, 
as seen from the numerical valu es of —Pm for t he first 7 longitudinal oscillation modes 
calculated by the Mathematica ( Wolfram 19961) and shown in table [T] it remains below 
its non-magnetic counterpart up to I = 5 modes. 

Note that the r.h.s of (|3.49p has a simple pole (z -1 ) singularity at z = 0, which is due to 
the neglected inertial term in (|3.41[) . As discussed above, this term becomes relevant for 
\z\ < Cm" 1 ^ 2 , where it cuts off the singularity at z~ x ~ Cm 1 ^ 2 . This cut-off integrated 
in (|3.50j) over \z\ < Cm^ 1 ^ 2 , where PJ n {z) ~ z for the odd modes, results in the damping 



rate correction 0(Cm 



-1/2*1 



which is a higher-order small quantity. 



3.5. Viscous energy dissipation 

Viscous damping rate can be found in an alternative much simpler way by considering 
the energy balance following from the dot product of (|3.1[) and v, which integrated over 
the drop volume yields 

\d t [ v 2 dV+ [ (V-n)vds = - [ (2Cae 2 + Cmj 2 ) dV, (3.52) 
2 Jv Js Jv 

where the first and second term on the l.h.s. stand for the time-variation of kinetic and 
surface energies, while the terms on the r.h.s. with the rate-of-strain tensor (e)ij = 
h( v i,j + v j,i) an d the dimensionless current density j — — + v x e account for the 
viscous and ohmic dissipations, respectively. As estimated above, viscosity gives rise to 
the tangential current density ~ Ha -1 in the Hartmann layer of the thickness ~ Ha -1 
that according to (|3.52D produces the ohmic dissipation ~ Cm/ Ha 3 ~ Ca/Ha which for 
Ha ^ 1 is negligible with respect to the viscous dissipation ~ Ca. Note that although 
the contribution of the Hartmann layer to the normal stress balance is important, its 
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m= 1 2 345 67 
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70 
11 


135 


14 


1232 
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220 


29 


1339 


7 


1350 


19 


4930 


11 




47 



27 



Table 2. The viscous damping rates —poi f° r the first 6 transversal oscillation modes. 



contribution to the energy dissipation is still negligible. This fact results in a substantial 
simplification of the solution procedure for the viscous damping rate. 

Thus, neglecting the ohmic dissipation and averaging the rest of (|3.52| over the period 
of oscillation and taking into account that the mean kinetic and surface energies for small 
amplitude harmonic oscillations are equal, we obtain a si mple expression for the v iscous 
damping rate in terms of inviscid leading order solution (jLandau fc Lifshitz 1119871 ) 



An =- / \io\ 2 dV/ 

Jv Jv 

For the longitudinal modes, this equation takes the form 

J [(rz^dzWof + (muiQ i 'r) 2 ]z 2 dz 



An 



(3.53) 



(3.54) 



fa wfz 2 dz 

Substituting tog(z) = ^B^P{ n (z)/z from (j3~2TJ)) into ([5751)> . after some algebra the last 
can be shown to be equivalent to (I3.50p . 

This approach is particularly useful for the transversal modes for which the conven- 
tional eigenvalue perturbation solution becomes excessively complicated and, thus, it is 
omitted here. In this case, using (|3.30[) we can represent (|3.53jl in terms of scalar potential 



0oi = - 

Substituting (p%(z) 
we obtain 



Jo (rdr(r 1 d r (pl)+m 2 (pl/r) 2 + {2mB r {^/r)y 



Jo [{dr^Y + (mpZ/r) 



'■ dz 



(3.55) 



l ^R e P{ n (z) from (pT28f into ([3755]) . after a lengthy algebra 



^ = -(2/+!) 



1(1 + l)(l - 2) - m 2 (l - 3) + (I 2 - m 2 )!^ 



(3.56) 



2(1(1 + 1) - m 2 ) 

where I™i is defined by ()3.51j) . Note that for 2D modes, defined by m = I, which 
are not affected by th e magnetic fi eld, we recover the well-known non-magnetic result 
/3§! = -(21 + l)(l - 1) (|Lamblll993h . For ot her indices (13351 can be verified by a direct 
integration of (|3.55p using the Mathematica ( Wolfranii~| [l9961 . As seen from the numerical 
values shown in table [5] the next even mode with m — I — 2 has the viscous damping 
rate which is by the factor of (I — 2)/(l — 4/5) lower than the non-magnetic counterpart 
given by m — I. Only for the modes with m < I — 4, the viscous damping rate in the 
magnetic field becomes higher than that without the field. 



The approach above is not directly applicable to the axisymmetric transversal modes 
which, as discussed at the end of §3.1.2[ are stationary in the leading-order inviscid ap- 
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proximation. For these overdamped modes, a flow with the velocity ~ 1/ Cm relative 
to the leading-order radius perturbation appears only in the first-order approximation, 
which according to (|3.52[) produces the same order ohmic dissipation. In this case, dissi- 
pation takes place on the account of the surface energy reduction, while that of the kinetic 
energy is negligible because it is by ~ 1/ Cm 2 smaller than the former. The contribution 
of the viscous dissipation in (|3.52[) is ~ Ca/ Cm 2 , which for a low viscosity and a high 
magnetic field is much smaller than the ohmic dissipation ~ 1/ Cm, and, thus negligible 
with respect to the latter. 



4. Conclusion 

In the present study, we have considered small-amplitude oscillations of a conducting 
liquid drop in a uniform DC magnetic field. Viscosity was assumed to be small but 
the magnetic field strong. Combining the regular and matched asymptotic expansion 
techniques we obtained a relatively simple solution to the associated eigenvalue problem. 
Firstly, we showed that the eigenmodes of shape oscillations are not affected by strong 
magnetic field. Namely, they remain the spherical harmonics as in the non-magnetic case. 
Strong magnetic field, however, constrains the liquid flow associated with the oscillations 
and, thus, reduces the oscillations frequency by increasing apparent inertia of the liquid. 
In such a field, liquid oscillates in a two-dimensional (2D) way as solid columns aligned 
with the field. Two types of oscillations are possible: longitudinal and transversal to the 
field. Such oscillations are weakly damped by strong magnetic field - the stronger the 
field, the weaker the damping, except for the axisymmetric transversal and 2D modes. 
The former are magnetically overdamped because the incompressibility constraint does 
not permit an axially uniform radial flow. The latter, which are transversal modes defined 
by the spherical harmonics with equal degree and order, I = m, are not affected by the 
magnetic field because these modes are naturally invariant along the field. In a uniform 
magnetic field, no electric current is induced and, thus, no electromagnetic force acts 
on such a 2D transversal flow because the associated e.m.f. is irrotational. Because the 
magnetic damping for all other modes decreases inversely with the square of the field 
strength, the viscous damping may become important in a sufficiently strong magnetic 
field. Consequently, the relaxation of axisymmetric transversal modes, whose viscous 
damping is negligible relative to the magnetic one, can be used to determine the electrical 
conductivity, while the damping of / = m modes can be used to determine the viscosity. 
The damping of all other modes is affected by both the viscous and ohmic dissipations. 
Although the latter reduces inversely with the square of the field strength while the 
former stays constant, an extremely strong magnetic field may be required for the viscous 
dissipation to be become dominant. 

As an example, let us consider a drop of Nickel of 1 cm in diameter (Rq = 5 x 10 -3 m) 
which, at the melting point (1455°C), has the surface tension 7 = 1.8 N/m, density 
p = 7.9 x 10 3 kg/in 3 , the dynamic vis cosity 77 = 4.9 x 10~ 3 N s/m 2 and the elec- 
trical conductivity a = 1.2 X 10 6 S/m (jGale k. Totemeier I [20041 ). The capillary time 
scale and frequency of non-magnetic fundamental mode (/ = 2) for such a drop are 
to = t/RJp/i ~ 23 ms and / = Jl(l - l)(l + 2)/(27rrn) w 19 Hz, respectively. The vis- 
cous damping time without the magnetic field ( Lamb II 1993^ 1 is t v /((21 + 1)(Z — 1)) m 8 s, 
where t v = pRq/t) w 40 s is the viscous time scale. Note that weak- viscosity approx- 
imation is applicable in this case because Ca = tq/t v = 5.8 x 10~ 4 is small. In the 
magnetic field of B = 5T, for which Cm = a B 2 Rq j \J p7 Rq « 87 > 1, the oscillation 
frequency of longitudinal fundamental mode (I, m) = (2,1) drops according to equation 
(j3T8| to f£ A = l)(/ + 2)/(27rr ) « 14 Hz, The corresponding viscous damping 
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time increases by the factor of two to — t v //3qi ~ 16 s, where — /3qi = 5/2 according to 
table Q] The magnetic damping time of this mode, for which (|3.32|) yields j3° — —4, is 
— ToCm/ ftl 0.5 s. According to this formula, for the magnetic damping time to exceed 
the viscous one, a magnetic field of B > 20 T is necessary. The relaxation time for the 
axisymmetric fundamental mode (I, m) — (2,0), which is magnetically over-damped, is 
—t Cm / (3f rj 28 ms, where j3f — 72 follows from (|3.37p . The magnetic field affects neither 
the frequency nor the damping rate of (I, m) — (2, 2) transversal oscillation mode, which 
is naturally invariant along the field. For the same reason, there is no magnetic damp- 
ing of this mode either. The first oscillatory transversal mode is (l,m) — (3,1) whose 
frequency drops according to (|3.27|) from /; = — 1)(Z + 2)/(27rr ) ~ 38 Hz without 

the magnetic field to = <J ly+i^J ^ I (27ttq) ~ 5 Hz in a strong magnetic field. The 
magnetic damping time for this mode in a 5T magnetic field is —tqCiti/ /3f rs 36ms, 
where j3f = 6800/121 follows from (|3.37p . The viscous damping time for this mode is 
— Tjz/ySoi ~ 6 s, where /3qi = 70/11 follows from table [5J The viscous damping is small 
relative to the magnetic one for this mode, and a magnetic field of about 65 T would be 
necessary for the magnetic damping time to become as long as the viscous one. 

In conclusion, this theoretical model provides a basis for the development of new mea- 
surement method of surface tension, viscosity and electrical conductivity of liquid metals 
using oscillating drop technique in a strong superimposed DC magnetic field. 

The author would like to thank Agris Gailitis and Raul Avalos-Zuniga for constructive 
comments and stimulating discussions. 
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